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Abstract 



A benchmark is set on the three-nucleon photodisintegration calculating the 
total cross section with modern realistic two- and three-nucleon forces (AV18, 
UrblX) using both the Faddeev equations and the Lorentz Integral Transform 
method. This test shows that the precision of three-body calculations involv- 
ing continuum states is considerably higher than experimental uncertainties. 
Effects due to retardations, higher multipoles, meson exchange currents and 
Coulomb force are studied. 
PACS: 21.45.+V, 25.20.Dc 
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I. INTRODUCTION 



The total photodisintegration cross sections of the three-body nuclei are important quan- 
tities for understanding electromagnetic processes in nuclei. Several approximate theoretical 
calculations have been performed in the past with simple NN-potential models The 
first calculations with complete final state interaction were carried out with the Lorentz In- 
tegral Transform (LIT) method using semirealistic || and realistic NN and 3N forces ||[7|]- 
In this paper we intend to set a benchmark on the results of the three-nucleon total photo- 
disintegration by calculating the cross section with modern realistic two- and three-nucleon 
forces. To this end we use two different methods, namely the Faddeev and the LIT ap- 
proaches. On the one hand the comparison between the results gives an idea of the precision 
reached by three-body calculations involving continuum states while on the other hand it 
may stimulate more accurate measurements. 

It has long been believed that the photonuclear cross section is dominated by unretarded 
dipole transitions. The contributions of magnetic and higher order electric multipoles, of 
retardation, as well as of meson exchange currents have been discussed extensively. However, 
rigorous results have only been given for the deuteron B-ID|. In this paper we give for the 
first time a quantitative discussion of these issues for the three-body nuclei. 

The paper is organized as follows. In Sec. 2 we briefly summarize the two methods for 
obtaining the total photodisintegration cross section. In Sec. 3 we present the comparison 
between the two calculations and discuss the effects of the Coulomb force, as well as the 
contributions beyond the unretarded dipole operator. Conclusions are drawn at the end. 



II. THEORETICAL FRAMEWORK 

The total unpolarized photodisintegration cross section is given by 

a T (E 7 ) = Rt(Ej) i-4vr 2 ^ . (1) 

In terms of the matrix elements of the transverse current operator J T the transverse response 
function Rt^E^ is 

R T (E,) = Jdf\^ f \J T ^ )\ 2 5(E,-E f + E -^-), (2) 

where \l/o is the three-nucleon bound state wave function with energy E$, and are 
final state wave functions. The term E^/(2M t ) in the energy conserving delta-function 
represents the small recoil energy of the whole 3-body nucleus with mass M t . In the following 
subsections it is described how Rt is calculated using the Faddeev equations and the LIT 
method. 



A. The Faddeev equations 



The way the transverse response function R? is calculated in the Faddeev framework for 
inclusive electron scattering without a 3N force is described in Refs. [[[]]] and [|12|]. Here we 
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sketch the corresponding derivation for the case when a 3N force is additionally taken into 
account. 

For a general operator O and using closure the function in Eq. (0) is rewritten as 



H 



(3) 



where E = Eq — E^/(2M t ) + E 1 and H is the 3N Hamiltonian. Next we define an auxiliary 



state 



1 



O|*o) 



E + ie-H 

which fulfils the following equation 

(E + ie-H) |*) = 0\%). 

This leads to an intermediate equation 

|*> = G (V 1 + V 2 + V 3 + V 4 )\V) + G O\^ ). 

Here Vi is the NN potential (we use the notation V\ = V23 etc.), V4 is a 3N force and Go the 
free 3N propagator. Both the operator O and the 3N force V4 can be written as a sum of 
three parts having the same symmetry under particle exchanges. 



(4) 

(5) 
(6) 



O 



1=1 
3 



(7) 



v 4 = J2v, 







i=l 



Introducing the Faddeev decomposition of the state |\&) 

I*) = GoEl^) (8) 

i=l 

we obtain from Eq. (|5|) (for three identical nucleons) the following equation on \Ui): 

(l-ViGo)^) = ViGoPlU,) + v} 1) G (l + P)\U 1 ) + Oi|* >, (9) 

where P is the sum of a cyclical and anticyclical permutation of 3 particles. Using the 
identities for the NN t-operator t\ 



(1-V 1 G )- 1 V 1 G = hG 
we obtain the final Faddeev-like equation on \Ui): 

\U X ) = (1 + ^0)0x1*0) + (txGoP + (l + tiG )U 4 (1) (l 
The response function is then obtained as 

R T (E,) = -- 9 (* \O\ (1 + P) G \U t ), 

71 



P)Go) |C/i). 



(10) 

(11) 
(12) 



where O is taken as the transverse current operator of Ref. [15] 

The response function Rt can also be obtained by means of direct integrations in Eq. (0). 
This gives not only a possibility for a check of numerics, but provides information about 
two- and three-body parts of the total cross section. 
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B. The Lorentz Integral Transform Method 



In the LIT approach the need to compute final state continuum wave functions is avoided 
l~3fl . In fact the response function Rt is obtained via evaluation and subsequent inversion 



of its LIT 

r°° R(E ) 
L(a R ,a!)= dE,— 1 ?!_ 2 . (13) 

JE th {Ery -CT R y+ af 

Neglecting the small recoil term in Eq. (0) E 7 = Ef — E and applying closure the transform 
L(cr) of the response Rt is found as 

L(vr, cri) = ($((Tr,(Ti)\*((Tr,(Ti)), (14) 

where \I/ is the localized solution of the Schrodinger-like equation 

(H - E - a R + iaj)^ = Q (15) 

with the source term Q = Jt^o- The wave function \l/ is the ground state solution for the 
same hamiltonian. 

In order to calculate \&o and ^ expansions on correlated hyperspherical harmonics (CHH) 
are used 

* = d^afa, (16) 

i 

where Q is a correlation operator and (pi is a totally antisymmetric basis set constructed 
from a spatial part Xi,n an d a spin-isospin part 9^ 

& = E ■ ( 17 ) 

The operator Q is a state dependent correlation operator. Further details can be found in 
Ref. 0. Calculating L(a R ,ai) for a sufficient number of a R and fixed 07, one obtains Rt 
from the inversion of the LIT (for details of the inversion see Ref. fT4|| ) 



C. Comparison between the two Approaches 

It is evident that the Faddeev and LIT methods are completely different. The Faddeev results 
are obtained essentially by solving Eq. ([11]) in momentum space while the LIT results are 
obtained by first solving Eq. (|15D in configuration space and then inverting Eq. (|l^) with 
L(a) given by Eq. flTjj). Therefore the comparison between the results obtained with the 
two methods is a very significant test of accuracy reached by these approaches. 
There is an interesting similarity between the two methods, namely, the fact that closure 
plays an important role in both approaches. Indeed it is just the use of closure that leads 
to Eqs. (^) and fli"5D, respectively. Actually these two equations are equal provided that 
<jj — > e. However in the LIT case the finite 07 in Eq. (|15|) effectively makes it a bound state 
problem, while the Faddeev-like equation (|TT|) remains a continuum problem. 
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In judging the quality of agreement between the results of the two methods one has to take 
into account that the LIT results are obtained in the "unretarded dipole approximation" . In 
the following we recall this approximation briefly. It consists in replacing the total current 
operator J(Q) by its limit J(Q = 0). Then one makes use of the continuity equation 



Q • J/o(Q) = £ 7 p/ (Q) 



(18) 



where p(Q) represents the charge operator. Expanding p/ (Q) in powers of Q and taking 
into account that p/o(<3 = 0) = 0, one gets 



Q • J f0 {Q = 0) = E 7 Q 
As the direction of Q is arbitrary, one has 



V Q p 



/o 



J /0 (Q = 0) = £ 7 [(V q p) 



fa 



Since the gradient of p for Q = is the dipole operator one writes 



/o 



fo 



(19) 



(20) 



(21) 



Thus Eq. ([I]) becomes the well known result for the photonuclear cross section in unretarded 
dipole approximation 



a = 4ii 2 —E~R D (E~ 

he 



with 



R D {E 1 )= Jdf\(* f \D\*o)\ 2 6(E 7 -E f + Eo 



(22) 



(23) 



For the LIT only R D is calculated. Since our energy range of interest is below pion threshold 
the unretarded dipole approximation should be very good. In the following discussion a check 
of this assumption will be presented in the framework of the momentum space Faddeev 
approach. In this framework one can use the representation of the electric multipoles of 
the transition operator obtained in Ref. |15|] employing the continuity equation (|Tj|) with 
Ery = Q: 



TfAQ) 



rpp , rpres 

1 JS, + 1 J£ ■> 



rpres 

1 Ji 



47T 



47T 



'J+ 1 



J 



dQYjs(Q) Pf . (Q) 



'2J + 1 



J 



dQ (Yij+iAQ) ■ J /o(Q) 



(24) 
(25) 

(26) 



Here ^ = ±1 and Yj J+11 are vector spherical harmonics. In the limit Q — > the contri- 
bution (p^) behaves as Q J+1 , while the contribution ( p5|) behaves as Q J . Expressing Tjg in 
terms of the density matrix elements pf o(R) one finds that 
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lim 



Q~ J Tjm(Q) 



lim 



Q~ J T p JM {Q) 



—i 



'J + i 



J (2J + 1)!! 



dRR J Y JM {R)p f0 (R) , 
(27) 



which is known as Siegert's theorem, namely that the electric multipole transitions are 
completely determined via the charge density in the limit Q — > 0. In particular, for J = 1 
one has 



(D 



M. 



fO 



6tc lim 



Q 



(28) 



This is the relation which has been used in the framework of the Faddeev calculation to test 
the quality of the unretarded El approximation. 

Corrections to the unretarded El approximation consist of two parts. The first part is 
due to the difference between T( M {Q) of Eq. ( |2"oD and its limit for Q — > 0. We shall refer to 
this correction as a retardation correction. The second part is Tj| s of Eq. (p6|). It is seen 
from above that both corrections vanish at small Q. Additional contributions come from 
magnetic multipoles (see Eq. (4) in ||15||) and higher EJ multipoles. 



III. RESULTS AND DISCUSSION 

For the calculation of the total photoabsorption cross section we use the AV18 NN 
interaction ||16|| , the 3-body force Urbana IX |17| and for the LIT case also the Coulomb 



force. Only the unretarded dipole approximation Ftp is considered within the LIT method. 
The Faddeev results also include retardation effects, additional multipoles of the usual one- 
body current and explicit MEC. 

We start the discussion with the benchmark test. In Fig. 1 the comparison between 
LIT and Faddeev results is presented in unretarded dipole approximation without (Fig. la) 
and with (Fig. lb) the 3N force. There are two curves for the LIT case showing the small 
uncertainties in the inversion of the transform. Besides small differences at higher energies 
the agreement between the two methods is excellent for AV18. Inclusion of the 3N force 
gives slightly poorer agreement in the peak region and at the low energy side of the peak, 
where the two curves have a relative shift of 0.2 MeV. We do not think that this small 
difference is of relevance at present; only in case of future high precision experiments would 
it be reconsidered. Fig. lb shows also the result for AV18 only. As already pointed out in 
Refs. H|7j one finds a decrease of the peak and enhancement of the tail due to the 3N force. 

In Fig. 2 we show the retardation effects (Fig. 2a) and those due to the additional 
contributions of T^ s and of magnetic and higher electric multipoles (Fig. 2b). Only the 
one-body current was employed to calculate these additional contributions. One sees that 
the retardation effects reduce the cross section. They are very small below 50 MeV (1 % or 
less) and become somewhat more sizable at higher energies (5 % at 100 MeV, 17 % at pion 
threshold). As seen in Fig. 2b the increase due to additional multipoles is rather similar with 
and without 3N force. Only beyond 80 MeV does the difference become somewhat more 
pronounced. At low energy there is a tiny decrease due to the additional El contribution. 
Above 30 MeV the higher multipole contributions are increasing up to 30 % at pion threshold. 
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It is interesting to note that one finds a partial compensation of retardation and higher 



multipole effects as predicted by Gerasimov via sum rule considerations [18[ . 

In Fig. 3 we illustrate the effect of the Coulomb force in the threshold region. The neglect 
of the Coulomb force leads to an increase of the cross section close to threshold (about 10 % 
at 7 MeV), but becomes rather small beyond 15 MeV. Considering the Coulomb force only 
in the ground state is not a good approximation. It is better to neglect the Coulomb force 
completely. 

Now we turn to the discussion of MEC effects. Generally it is a problem to construct 
a consistent exchange current for a given NN potential. Using, however, minimal coupling 
one obtains a MEC with the correct divergence of the current also in case of realistic NN 
potentials P,|l9"ll, while the rotor of the current remains model dependent. Here we use only 
the standard n- and p-like exchange currents for AV18, which are determined according to 



the Riska prescription |fL9|| . In the present Faddeev calculation MEC are added to the one- 
body current without making a multipole decomposition, therefore without using Eq. (f24"|). 
This differs from the usual approach where the gauge condition Eq. is applied to obtain 
a representation similar to that of Eq. (^) and thus only MEC contributions beyond it are 
calculated explicitly. Such a procedure has the advantage that one violates gauge invariance 
in the minimal way since the divergence of the full current is taken into account correctly. 
For example, in a calculation with 3N forces one automatically takes into account a large 
part of three-nucleon MEC effects. Fig. 4 illustrates that the result with explicit MEC 
differs rather substantially from that with one-body current and implicit MEC via Siegert's 
theorem. This suggests that the MEC are not fully consistent with the potential. Indeed, 
using Eq. for a current satisfying the continuity equation the explicit MEC contribution 



to the El multipole would manifest itself only in the residual term (26) which vanishes for 
Q — > 0. Even away from this limit TJ| S is probably small, as it is in the case of the deuteron 
|20| . In addition, MEC contributions to other multipoles could not decrease the cross section, 
since the cross section is the incoherent sum of all multipoles and without explicit MEC such 
multipole contributions are negligible at low energies (see Fig. 2b). 

The contributions of the two- and three-body break up channels are shown in Fig. 5a. 
The three-body break up cross section becomes larger than that of the two-body break 
up already at 14 MeV and at higher energies it is the dominating channel. In Fig. 5b we 
compare the total three-body break up cross section with that of the final isospin T=3/2 
channel (which is three-body break up exclusively). It is interesting to see that at low energy 
the T=l/2 channel gives a small contribution to the three body break up, while it becomes 
increasingly important at higher energies. 

In Fig. 6 we compare the triton results with experimental data. One finds a rather 
good agreement between theory and experiment at low energy for both channels. Since the 
experimental situation at higher energy is not settled, a definite comparison between theory 
and experiment cannot be made. In fact the two-body break up data are rather scattered 
and there is limited experimental information on the three-body break up. 

In Fig. 7 we compare the total 3 He cross section with data (the results for 3 H with AV18 
and UrblX are presented in 0). One finds a rather good agreement between theoretical 
and experimental data and it is evident that the 3N force improves the comparison with 
experiments at low energy. On the other hand due to the insufficient precision of the 
experimental data a more conclusive comparison between theory and experiment cannot be 
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made. 

We summarize our results as follows. Benchmark results of high precision have been 
set for the three-nucleon total photoabsorption cross sections. The theoretical quality of 
the results is considerably higher than the experimental uncertainties. Of course, one has 
to take into account that these experimental data are rather old. Thus new experimental 
activity in this field would appear to be timely. The classical approximation for photonuclear 
cross sections, i.e. the unretarded dipole approximation, is extremely good below 50 MeV. 
It is also rather good at even higher energies (error at 100 MeV below 20 %), since one 
has a partial cancellation of the El retardation effects and contributions from additional 
multipoles. The Coulomb force effects on the cross section become negligible about 6-7 
MeV above threshold but are important at lower energy. 
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FIGURES 



FIG. 1. Comparison of the Faddeev and LIT results for the total 3 H photoabsorption cross 
section in unretarded dipole approximation without (a) and with (b) 3N force. The dots are the 
Faddeev results and the two curves represent the bounds for the inversion of the LIT. The dotted 
curve in (b) is the result with AV18 only. 

FIG. 2. (a) Retardation effects to the unretarded El triton cross section [p re % — a unre t) / 'cr U nret 
and (b) retardation and additional contributions (see text) (a to t — cr unre t)/a unre t. 

FIG. 3. Relative effect of the Coulomb force on the 3 He photoabsorption cross section: no 
Coulomb at all, i.e. triton result with shifted threshold (dashed curve), no Coulomb in final state 
interaction (dotted curve). 

FIG. 4. 3 He photoabsorption cross section without a 3N force for the full one-body current and 
implicit MEC contribution via application of Eq. (24) (magnetic one-body transitions considered 
via Eq. (4) of Ref. [15]) (full dots) and result for the full one-body current with additional explicit 
MEC but without application of Eq. (24) (empty dots). Curves are spline interpolations. 

FIG. 5. (a) Two-body (dotted curve), three-body (dashed curve) and total (full curve) break up 
cross sections of the 3 H photodisintegration (full one-body current and implicit MEC contribution 
as in Fig. 4). (b) Comparison of the total three-body break up cross section (dashed curve) with 
the cross section in the final state isospin channel T=3/2 in the unretarded dipole limit. 

FIG. 6. Two-body (a) and three-body (b) break up cross sections for the 3 H photodisintegra- 
tion (full one-body current and implicit MEC contributions as in Fig. 4). The full curves represent 
Faddeev results. Experimental data in (a) from [21] (circles), [22] (squares), [23] (diamonds), [24] 
(triangles up), [25] (triangles down), [26] (X). The dotted lines in (b) are bounds of experimental 
data from [23]. 

FIG. 7. Total 3 He photoabsorption cross section. LIT results in unretarded dipole approxima- 
tion with AV18 alone (dash-dotted curve) and with AV18+UrbIX (full curve); the dotted lines are 
bounds of experimental data from [23] and the dots are data from [27]. 
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